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INTRODUCTION 

The  Department  of  Defense  (DOD)  has  numerous  facilities  engaged  in 
the  production  of  various  types  of  explosives  and  munitions  used  by  mil- 
itary services.  In  most  cases  the  production  of  ammunition  utilizes 
assembly  line  procedures.  Projectiles  pass  through  various  stages  of 
preparation:  filling  with  explosive,  fuzing,  marking,  and  packing. 

Hazardous  operations,  such  as  the  filling  of  the  projectile  case  with  an 
explosive  in  a powder  form  and  the  compaction  of  the  powder  by  hydraulic 
press,  are  accomplished  in  protective  cells  that  are  intended  to  confine 
the  effects  of  an  accidental  explosion. 

Most  of  the  existing  production  facilities  were  built  in  the  1940s. 
With  few  exceptions,  the  manufacturing  technology  and  existing  equipment 
represent  the  state  of  the  art  as  of  1940.  The  production  equipment  was 
operated  extensively  during  World  War  II,  again  during  the  Korean  con- 
flict, and  recently  during  the  Southeast  Asia  war.  Much  of  this  equip- 
ment and  the  housing  structures  have  been  operating  beyond  their  designed 
capacities  [1], 

DOD  is  conducting  an  ammunition  plant  modernization  program  [2]  that 
is  intended  to  greatly  enhance  safety  in  the  production  plants  by  pro- 
tective construction,  automated  processing,  and  reduction  of  personnel 
involved  in  hazardous  operations. 

In  1969  a joint-service  manual  [3]  was  published  to  provide  guidance 
to  the  structural  designers  of  munition  plants.  The  objectives  of  the 
manual  were  to  establish  design  procedures  and  construction  techniques 
to  prevent  propagation  of  explosions  from  one  building,  or  part  of  a 
building,  to  another;  to  prevent  mass  detonations;  and  to  provide  protec- 
tion for  personnel  and  equipment.  The  manual  establishes  blast-load 
parameters  for  designing  protective  structures,  provides  methods  for  cal- 
culating the  dynamic  response  of  concrete  walls,  and  establishes  con- 
struction details  for  developing  required  strength.  The  design  method 
accounts  for  close-in  effects  of  a detonation  with  its  associated  high 
pressures  and  nonuniformity  of  loading  on  protective  barriers.  A detail- 
ed method  for  assessing  the  degree  of  protection  afforded  by  a protective 
facility  did  not  exist  prior  to  this  manual's  publication;  consequently, 
the  manual  represents  a significant  improvement  in  design  methods.  The 
simplifications  made  in  the  development  of  the  design  procedures  have 
been  presented  in  the  manual.  The  analysis  of  a structure  using  the 
design  procedure  will  generally  result  in  a conservative  estimate  of  the 
structure's  capacity;  therefore,  structures  designed  using  these  proce- 
dures will  generally  be  adequate  for  blast  loads  exceeding  the  assumed 
load  conditions  [3], 
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Even  with  the  simplifications  presented  in  Reference  3,  the  compu- 
tational procedures  are  complex  and  time-consuming.  An  automated  pro- 
cedure was  required  to  give  structural  designers  the  capability  to 
perform  rapid  analysis  of  the  structural  safety  of  blast-resistant  con- 
struction. The  design  parameters  interact  in  a complex  way  since  the 
procedure  is  both  nonlinear  and  dynamic.  From  a design  point  of  view 
an  optimization  procedure  was  required  to  minimize  cost  and  maximize 
safety  since  blast-resistant  construction  has  been  reported  to  cost 
3 to  5 times  as  much  as  conventional  construction.  Therefore,  the  first 
objective  was  to  automate  the  analysis  procedures  for  determining 
structural  response  of  reinforced  concrete  slabs  having  a bilinear 
stiffness  representation  and  subjected  to  blast  shock  and  gas  pressures. 
Concrete  slabs  are  the  basic  element  forming  sidewalls,  roofs  and  floors 
of  cells  designed  to  confine  the  effects  of  accidental  explosions.  The 
second  objective  was  to  provide  an  optimum  design  procedure  for  laced 
and  unlaced  reinforced  concrete  slabs  that  will  automatically  produce 
a least-cost  design  for  a given  slab  geometry,  material  properties,  and 
explosive  weight  for  both  feasible  and  nonfeasible  starting  points. 


COMPUTER  PROGRAM 

The  computer  program  was  written  in  FORTRAN  IV  for  use  with  Control 
Data  6600  series  computers.  The  program  is  composed  of  four  areas: 

1.  Blast  Load  Determination 

2.  Structural  Analysis  Parameters 

3.  Dynamic  Response 

4.  Optimization 

The  blast-load  determination  is  accomplished  by  subroutines  BLA, 

PIC,  SGRID,  HBA,  RATIO,  GRID,  GAS  INTERP,  EQUIV,  HEDATA,  ARDC,  SHOCK, 
and  TNT.  The  subroutines  read  the  explosive  weight  and  type  and  cell 
geometry,  and  then  compute  the  equivalent  spherical  weight  of  TNT  and 
the  equivalent  pressure  loading  using  the  geometry  of  the  wall  and 
charge  location.  Both  the  shock  pressure  and  its  duration  and  the  gas 
pressure  and  its  duration  are  calculated.  Using  the  duration  and 
pressure  data  for  both  shock  and  gas,  the  program  computes  an  equivalent 
triangular  pressure  loading  for  each  part  and  adds  both  together  to 
produce  the  resultant  shown  in  Figure  1.  The  total  impulse  is  then 
determined. 

The  structural  analysis  is  accomplished  by  subroutines  SSTIFF,  LACE, 
DOOR  1,  DOOR  2,  DOOR  3,  DOOR  4 and  DOOR  5.  These  routines  compute  the 
stiffness,  resistance,  and  equivalent  mass  of  the  slab  using  input 
material  properties.  Both  flexure  and  shear  are  considered.  Openings 
(doors  or  windows)  in  walls  are  allowed. 
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The  dynamic  response  calculation  is  accomplished  in  subroutine  RESP. 


The  program  determines  the  response  of  the  slab  modeled  as  an  equivalent 
dynamic  single-degree-of-f reedom  system  with  bilinear  stiffness  and 
pressure  loading  as  shown  in  Figure  1.  The  solution  technique  is  based 
on  a Newmark  iteration  method. 

When  a thickness  of  sand  is  specified  for  composite  construction 
(i.e.,  two  slabs  with  sandfill),  the  program  computes  the  impulse  capac- 
ity of  the  first  slab  using  half  the  mass  of  the  sand  as  acting  with  the 
wall.  Figures  6-38  and  6-39  of  Reference  3 give  the  attenuation  of  the 
blast  wave  on  sand  for  evaluation  of  the  impulse  capacity  of  the  second 
wall.  The  optimization  of  an  initial  design  is  accomplished  in  sub- 
routines OPT,  MIN1MZ,  PMINZ,  DMINZ,  GETE,  SUMRY,  TLEFT , AND  GCOMP.  The 
methodology  used  is  that  of  a penalty  function  with  individual  mini- 
mization sequences  being  accomplished  by  the  Powell  method. 

Program  Input 


The  program  input  consists  of  five  or  six  cards  per  case.  Additional 
cases  can  be  grouped  together.  Two  blank  cards  are  used  after  the  last 
case.  The  users  guide  contained  in  the  program  is  given  here  to  assist 
in  understanding  the  input.  Card  format  is  8F10.0  except  as  noted. 
Figure  2a  is  an  input  data  sheet  to  be  used  in  conjunction  with  Figures 
2b  and  2c,  which  show  the  slab  geometry  and  orientation  that  must  be 
followed.  The  input  required  for  each  card  is  described  below. 


CARD  1 

COL  2 COL  68  HEADING 

COL  69  COL  70  OPTIMIZATION  = 1. 

COL  71  COL  72  FLAG  KQ  0 FOR  PRESSURE  CALCULATION  EQ  1. 

FOR  INPUT  PRESSURE 

COL  73  COL  74  FLAG  FOR  REBAR  DIAMETER  RATHER  THAN  AREA  = 1 

COL  75  COL  76  FLAG  FOR  IMPULSE  GRID  = 1 

COL  77  COL  78  DOOR  FLAG  = 1 

COL  79  COL  80  DOOR  EQUILIBRIUM  ITERATION  = 1 

OTHERWISE  = 0 

CARD  2 

COL  1 COL  10  WEIGHT  OF  ACTUAL  EXPLOSIVE  LB 

COL  11  COL  20  EXPLOSIVE  NUMBER  SEE  TABLE  2 

COL  21  COL  30  EXPLOSIVE  LENGTH/DIAMETER  RATIO 

COL  31  COL  40  PROJECTILE  CASE  WEIGHT/EXPLOSIVE  WEIGHT  RATIO 

COL  41  COL  50  AMBIENT  PRESSURE  PS1A  (DEFAULT  14.69  PSI) 

COL  51  COL  60  AMBIENT  TEMPER ATURE°C  (DEFAULT  20°) 

COL  61  COL  70  ALTITUDE  KFT  (WHEN  PRESSURE  AND  TEMPERATURE 

NOT  SPECIFIED) 

COL  71  COL  80  EFFECTIVE  IMPULSE  FRACTION  COMPOSITE  CONSTRUCTION 


CARD  3 
COL  1 

COL  10 

RA  DISTANCE  CHARGE  TO  WALL  FT 
OR  EQUAL  IMPULSE  PSI-MS  IF  FLAG  = 1.0 

COL  11 

COL  20 

H WALL  HEIGHT  FT 

COL  21 

COL  30 

EL  WALL  LENGTH  FT 

COL  31 

COL  40 

HL1T  HEIGHT  CHARGE  FT 

OR  EQUAL  PRESSURE  PS1  IF  FLAG  = 1.0 

COL  41 

COL  50 

ELLIT  DISTANCE  CHARGE  TO 
LEFT  SIDE  WALL  FT 

OR  EQUAL  DURATION  (MS)  IF  FLAG  = 1.0 

COL  51 

COL  60 

CELL  VOLUME  FOR  GAS  PRESSURE 

COL  61 

COL  70 

CELL  VENT  AREA  FOR  GAS  PRESSURE 

COL  71 

EQ1  FOR  FLOOR  REFLECTION 

COL  72 

EQ1  FOR  ROOF  REFLECTION 

COL  73 

EQ1  FOR  LEFT  WALL  REFLECTION 

COL  7 4 

EQ1  FOR  RIGHT  WALL  REFLECTION 
OTHERWISE  EQ  0 

CARD  4 
COL  1 

COL  10 

FC  DYNAMIC  CONCRETE  STRESS  PSI 

COL  11 

COL  20 

FY  DYNAMIC  STEEL  STRESS  PSI 

COL  21 

CCL  30 

TC  THICKNESS  CONCRETE  IN. 

COL  31 

COL  40 

THETA  ALLOWABLE  ROTATION  DEGREES 

COL  41 

COL  50 

NSIDE  NUMBER  OF  SIDES  WALL  FIXED 

1.0 

BOTTOM  SIDE  FIXED 

2.0 

BOTTOM  AND  SIDE  FIXED 

3.0 

2 SIDES  AND  BOTTOM  FIXED 

4.0 

4 SIDES  FIXED 

5.0 

SIMPLE  SUPPORTED  BEAM  FIXED  AT  TOP  AND  BOTTOM 

6.0 

FIXED  BEAM  AT  TOP  AND  BOTTOM 

7.0 

BEAM  BOTTOM  FIXED  TOP  SIMPLE 

COL  51 

COL  60 

TSAND  SAND  THICKNESS  FT 

COL  61 

COL  70 

Bl.  LACING  LENGTH 

COL  71 

COL  80 

SL  LACING  SPACING  IN. 

CARD  5 
OPTION  2=0 


COL  1 

COL  10 

ASVT  AREA  VERTICAL  STEEL  BLAST  SIDE/FT 

COL  11 

COL  20 

ASVB  AREA  VERTICAL  STEEL  OPPOSITE  SIDE/FT 

COL  21 

COL  30 

ASHT  AREA  HORIZONTAL  STEEL  BLAST  SIDE/FT 

COL  31 

COL  40 

ASIIB  AREA  HORIZONTAL  STEEL  OPPOSITE  SIDE/FT 

COL  41 

COL  50 

DVT  DEPTH  TO  VERTICAL  STEEL  BLAST  SIDE  IN. 

COL  51 

COL  60 

DVB  DEPTH  TO  VERTICAL  STEEL  OPPOSITE  SIDE 

IN. 

COL  61 

COL  70 

DHT  DEPTH  TO  HORIZONTAL  STEEL  BLAST  SIDE 

IN. 

COL  71 

COL  80 

DIIB  DEPTH  TO  HORIZONTAL  STEEL  OPPOSITE  SIDE  IN 

OPTION  2 = 1 

COL  1 

COL  10 

BAR  SIZE  VERT  OPPOSITE  SIDE 

COL  11 

COL  20 

BAR  SIZE  VERT  BLAST  SIDE 

COL  21 

COL  30 

BAR  SIZE  IIORIZ  OPPOSITE  SIDE 

COL  31 

COL  40 

BAR  SIZE  IIORIZ  BLAST  SIDE 
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COL  41 

COL  50 

BAR  SPACING  VERT  OPPOSITE  SIDE 

COL  51 

COL  60 

BAR  SPACING  VERT  BLAST  SIDE 

COL  61 

COL  70 

BAR  SPACING  HOR1Z  OPPOSITE  SIDE 

COL  71 

COL  80 

BAR  SPACING  HORIZ  BLAST  SIDE 

CARD  6 

COL  1 

COL  10 

BAR  DEPTH  VERT  OPPOSITE  SIDE 

COL  11 

COL  20 

BAR  DEPTH  VERT  BLAST  SIDE 

COL  21 

COL  30 

BAR  DEPTH  HORIZ  OPPOSITE  SIDE 

COL  31 

COL  40 

BAR  DEPTH  HORIZ  BLAST  SIDE 

OPTION  4 = 1 
CARD  6 
BLAST  DOOR 

PARAMETERS 

DEPTH  FROM  O UTTER  CONCRETE  SURFACE 
TO  CENTER  OF  BAR 

COL  1 

COL  10 

DOOR  HEIGHT 

COL  11 

COL  20 

DOOR  WIDTH 

COL  21 

COL  30 

DISTANCE  FROM  LEFT  SIDE  TO  DOOR 

COL  31 

COL  40 

DOOR  REACTION  OR 

COL  41 

COL  50 

DOOR  RESISTANCE  FOR  CALCULATION  OF  REACTION 

COL  51 

COL  60 

DISTANCE  TO  FLOOR 

Note:  All  values  are  fixed  point,  except  for  reflection  code  and  options. 


The  explosive  number  (Card  2)  refers  to  the  list  of  explosives  in 
able  1.  This  is  used  to  compute  explosive  equivalence.  The  length/ 
diameter  ratio  for  an  explosive  sphere  is  0.0,  which  gives  a shape 
factor  of  1.0.  For  an  uncased  explosive  the  case  explosive  weight  ratio 
is  0.  For  sea  level  calculations  the  ambient  air  pressure,  Pamb’  and 
temperature,  Tam^,  and  altitude  can  be  left  blank  and  will  default  to 
14.69  psi  and  20°C.  If  the  flag  in  the  heading  card  is  set  to  1,  the 
impulse,  duration,  and  pressure  will  be  read  on  Card  3.  If  the  flag  is 
left  blank,  the  charge  to  wall  distance,  charge  height,  and  distance 
from  the  left  side  will  be  read.  If  NSIDE  is  left  blank,  the  program 
will  sum  the  number  of  reflecting  sidewall  surfaces  specified  on  Card  3. 
The  separate  use  of  NSIDE  is  helpful  when  a frangible  wall  is  present, 
which  creates  a shock  reflection  but  does  not  provide  any  support. 

When  optimization  and  composite  construction  are  specified  together, 
the  program  will  optimize  the  design  to  resist  the  given  or  computed 
impulse.  For  the  case  when  two  walls  are  acting  together  each  resisting 
a portion  of  the  impulse  it  is  necessary  to  specify  the  effective  impulse 
to  be  applied  to  the  wall  under  design.  The  total  impulse  is  multiplied 
by  the  decimal  number  specified  on  Card  2. 
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Table  1.  List  of  Explosives 


Explosive 

Number 


Explosive  Name  and  Composition 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


TNT 

TNETB 

EXPLOSIVE  D 

PENTOLITE  (PETN/TNT  50/50) 
PICRATOL  (EXPLOSIVE  D/TNT  52/48) 
CYCLOTOL  (RDX/TNT  70/30) 

COMP  B (RDX/TNT/WAX  59.4/39.6/1.0) 
RDX/WAX  (98/^) 

COMP  A- 3 (RDX/WAX  91/9) 
TNETB/AL  (90/10) 

TNETB /AL  (78/22) 

TNETB/AL  (72/28) 

TNETB/AL  (65/34) 

TRITONAL  (TNT/AL80/70) 
RDX/AL/WAX  (88/10/2) 
RDX/AL/WAX  (89/20/2) 
RDX/AL/WAX  (74/21/5) 
RDX/AL/WAX  (74/22/4) 
RDX/AL/WAX  (62/33/5) 

TORPEX  II  (RDX/TNT/AL  42/40/18) 

H6  ( RDX / TNT / AL / WAX  45/29/21/5) 
HBX-1  (RDX/TNT/AL/ WAX  40/38/16/5) 
HBX-3  (RDX/TNT/AL /WAX  31/29/35/5) 
TNETB/RDX/AL  39/26/35) 
ALUMINUM 
WAX 
RDX 
PETN 
TETRYL 
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Example  Problems 


The  first  example  is  a sidewall  of  a blast  cell  with  a roof.  The 
concrete  wall  is  32  feet  long,  12  feet  high,  and  2 feet  thick  with  4 feet 
of  sand  in  composite  action.  Note  that  half  of  the  input  thickness  of 
the  sand  is  used  by  the  program  as  added  mass  to  the  wall.  The  wall  is 
restrained  at  the  floor,  roof,  and  left  side;  the  right  side  is  free. 
Since  the  three-side-fixed  option  condition  assumes  the  sides  and  the 
bottom  to  be  fixed,  the  wall  must  be  reoriented  when  filling  out  the 
input  form  (Figure  3).  Thus,  a height  of  32  feet  and  a length  of  12 
feet  are  used  to  properly  orient  the  free  edge.  An  allowable  support 
rotation  of  12  degrees  is  given  which  assumes  lacing  reinforcement  will 
be  used. 

Figure  4 gives  the  results  of  the  analysis.  A blast  impulse  of 
2,406  psi-ms  was  determined.  The  section  properties  are  given.  Since 
shear  exceeds  the  allowable,  lacing  must  be  provided  for  the  difference. 
The  yield-line  location  is  given.  An  ultimate  resistance  of  101.9  psi 
and  a stiffness  of  825  psi  were  determined.  The  impulse  capacity  of  the 
wall  is  4,365  psi-ms,  which  is  much  greater  than  the  loading  of  2,406 
psi-ms;  this  indicates  the  design  is  conservative.  If  a second  wall  of 
the  same  construction  were  present  and  acted  with  the  first  in  composite 
construction.  Figures  6-38  and  6-39  of  Reference  3 could  be  used  to 
determine  its  impulse  capacity  and  the  total  capacity  of  both  walls; 
the  scaled  values  of  impulse,  sand,  and  concrete  thicknesses  are  used. 

Figure  5 gives  the  input  data  for  a roof  of  a blast  cell  32  by  15 
feet.  The  32-foot  side  is  used  as  the  height  to  agree  with  the  fixity 
condition.  Figure  6 presents  the  computer  analysis.  In  this  case  sand 
fill  is  not  present,  and  the  wall  response  is  calculated. 

A maximum  deflection  of  18.98  inches  was  determined;  this  can  be 
compared  with  the  allowable  12-degree-rotation  deflection  of  18.6  inches. 
Since  the  maximum  deflection  exceeds  the  12-degree-rotation  deflection, 
collapse  of  the  wall  is  indicated.  Average  and  maximum  scale  velocity 
are  given.  The  appendix  gives  two  additional  examples,  wherein  hand 
calculations  are  compared  with  computer  results. 

THEORETICAL  DEVELOPMENT 

Blast  Loads  and  Structural  Response 

In  general,  the  methods  used  in  the  computer  program  follow 
Reference  3,  and,  as  such,  the  accuracy  of  both  is  the  same.  Since 
these  are  discussed  in  detail  in  References  3 and  4,  they  will  not  be 
presented  here.  The  solution  of  the  dynamic  response  equation  of  motion 
has  been  found  to  agree  very  closely  with  the  response  chart  of 
Reference  3.  Additionally,  the  solution  covers  a wider  range  and,  thus, 
is  more  accurate  in  the  areas  not  defined  by  the  response  chart.  When 


rt 
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the  loading  is  less  than  one  hundredth  of  the  natural  period,  the 
response  is  determined  by  impulse  equilibrium.  The  basic  dynamic  model 
is  limited  to  one  mode  of  response  and  does  not  consider  higher  modes. 

The  ultimate  moment  capacity,  Mu,  of  the  slab  is  based  on  Equation 
5-4  of  Reference  3,  as  follows: 


M 

u 


A’  f 

+ — 5 — — (d-d') 
b 


where  = area  of  compression  reinforcement 

A = area  of  tension  reinforcement 
s 

b = width 

a = depth  of  equivalent  rectangular  stress  block 

f = design  steel  stress 

d = distance  from  extreme  compression  fiber 
to  centroid  of  tension  reinforcement 

d'  = distance  from  extreme  compression  fiber 
centroid  to  compression  fiber 

This  equation  for  equal  reinforcement  in  tension  and  compression 
reduces  to 


A'  f 

M = — 5 — - (d  - d') 
u 


The  action  of  the  concrete  in  compression  is  neglected,  because 
crushing  at  high  rotations  is  assumed  to  occur.  This  results  in  dis- 
engagement of  the  concrete  cover.  When  support  rotations  are  restricted 
by  lack  of  lacing,  this  equation  becomes  conservative.  However,  the 
more  conventional  concrete  analysis  procedures  were  not  included  to 
conform  with  the  methodology  given  in  Reference  3. 

The  blast  impulse  computation  is  restricted  to  a geometry  in  which 
the  slab  height-to-length  ratio  is  greater  than  0.2.  The  modification 
made  by  the  Naval  Surface  Weapons  Center  to  the  original  Picatinny 
Arsenal  Program  did  not  affect  the  results  significantly  for  most  cases. 
However,  it  did  remove  several  minor  problem  areas,  such  as  the  location 
of  the  charge.  The  blast  impulse  has  all  the  limitations  associated 
with  the  original  Picatinny  programs  that  are  caused  by  limitations  in 
the  test  data.  It  assumes  the  charge  is  an  equivalent  sphere  of  TNT. 
Shape  effects,  explosive  equivalence,  and  explosive  casings  are  con- 
sidered, but  only  in  an  empirical  manner  as  a result  of  limited 
available  data. 
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Structural  Optimization 


The  optimization  problem  consists  of  finding  the  least-cost  struc- 
ture that  satisfies  all  the  design  constraints;  or,  stated  in  opti- 
mization terms:  ^ 

Find  such  that  M(X)  is  a minimum  and 

g±(X)  <0  i = 1,  2,  N 

where  X = vector  of  design  variables 

N = number  of  design  constraints 

g = vector  of  design  constraints 

M = objective  function 

Specifically  for  this  problem,  the  design  variables  selected  are  areas 
of  steel  reinforcement  and  thickness  of  concrete.  The  design  constraints 
are  the  flexural  and  shear  limits.  The  objective  function  consists  of 
the  costs  of  formwork  and  concrete  flexural  and  shear  reinforcement. 

Fixed  Variables 


W 

H 

EL 

h 

l 

R 

a 

I 


0 


explosive  weight 
wall  height 
wall  length 

height  of  explosive  above  floor 
distance  of  explosive  from  left  side  of  wall 
distance  of  explosive  from  wall 
reflection  code 

ultimate  dynamic  concrete  strength 
dynamic  yield  strength  of  reinforcing  steel 
rotation  criterion 


Design  Parameters,  X 


t concrete  thickness 

c 

X = AV  area  of  vertical  reinforcing  steel 

AH  area  of  horizontal  reinforcing  steel 

Constraints,  g (X) 


6(X)  = 6(0),  maximum  deflection 

V(X)  < VC  for  0 < 2 deg,  maximum  shear 


t >12,  minimum  thickness 

c “ 

AV  > 0.0025  bd 

AH  > 0.0025  bd 


minimum  steel  reinforcement 


The  methodology  [5,6]  selected  uses  the  unconstrained  minimization 
approach.  The  problem  is  converted  to  an  unconstrained  minimization  by 
constructing  a function,  <(>,  of  the  general  form 


<KX,  r)  = M(X)  + P[gl(X), 


...,  gn(X),  r] 


For  this  problem  the  interior  penalty  function  technique  was  selected. 
This  methodology  is  suitable  when  gradients  are  not  available,  and, 
because  the  method  uses  the  feasible  region,  a useable  solution  always 
results.  The  objective  function  is  augmented  with  a penalty  term  that 
is  small  at  points  away  from  the  constraints  in  the  feasible  region, 
but  increases  rapidly  as  the  contraints  are  approached.  The  form  is  as 
follows: 


4>  (X,  r)  = M(X)  - r 


N 

Z 

3=1 


gjOO 


where  M is  to  be  minimized  over  all  X satisfying  g (X)  < 0,  j = 2 ...  N. 
Note  that  if  r is  postive,  then,  since  at  any  interior  point  all  of  the 
terms^in  the  sum  are  negative,  the  effect  is  to  add  a positive  penalty 
to  M(X) . As  the  boundary  is  approached,  some  g (X)  will  approach  zero, 
and  the  penalty  will  increase  rapidly.  The  parameter,  r,  will  be  made 
successively  smaller  in  order  to  obtain  the  constrained  minimum  of  M. 


Objective  Function,  F 


Cost  = F = H • EL  • t 


+ (AV  + AH) (EL  • H)C  + 


(As) (EL 


H)Ct 


where  C = cost  of  concrete  ($/cu  ft) 
c 

C = cost  of  horizontal  and  vertical  reinforcement  ($/cu  in.) 
s 

C = cost  of  lacing  reinforcement  ($/cu  in.) 

L 

A = area  lacing  reinforcement  ($/cu  in.) 


= F + 


N 

r Z 
3=1 


8j<x) 


where 


penalty  parameter. 
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The  program  requires  a starting  point  in  the  feasible  region 
before  optimization  can  proceed.  This  is  accomplished  automatically 
by  the  program  by  incrementing  the  design  variables  until  a feasible 
point  is  reached. 

An  algorithm  which  comprises  the  steps  most  commonly  used  is  as 
follows : 

1.  Given  a starting  point,  Xq,  satisfying  all  g j (X)  < 0 and  an 

initial  value  for  r,  minimize  4>  to  obtain  X , . 

min 

2.  Check  for  convergence  of  Xmin  to  the  optimum. 

3.  If  the  convergence  criterion  is  not  satisfied,  reduce  r by 
r ■*-  rc,  where  c < 1. 

4.  Compute  a new  starting  point  for  the  minimization,  initialize 
the  minimization  algorithm,  and  repeat  from  step  1. 

The  logic  diagram  for  the  interior  penalty  functions  technique  is 
shown  in  Figure  7. 

The  minimization  for  tf>(X,  r)  shown  in  Figure  7 is  accomplished  by 
a method  developed  by  Powell  using  conjugate  directions  [5,6]. 

Powell's  method  can  be  understood  as  follows:  Given  that  the 

function  has  been  minimized  once  in  each  of  the  coordinate  directions 
and  then  in  the  associated  pattern  direction.  Discard  one  of  the 
coordinate  directions  in  favor  of  the  pattern  direction  for  inclusion 
in  the  next  m minimizations,  since  this  is  likely  to  be  a better 
direction  than  the  discarded  coordinate  direction.  After  the  next 
cycle  of  minimizations,  generate  a new  pattern  direction,  and  again 
replace  one  of  the  coordinate  directions.  This  process  is  illustrated 
in  Figure  8. 

Figure  9 is  a logic  diagram  for  the  unconstrained  minimization 
algorithm.  The  pattern  move  is  constructed  in  block  A,  then  used  for 
a minimization  step  (blocks  B and  C) , and  then  stored  in  Sn  (block  D) 
as  all  of  the  uirections  are  up-numbered  and  is  discarded.  The 
direction  Sn  will  then  be  used  for  a minimizing  step  just  before  the 
construction  of  the  next  pattern  direction.  Consequently,  in  the 
second  cycle,  both  X and  Y in  block  A are  points  that  are  minima  along 
Sn,  the  last  pattern  direction.  This  sequence  will  impart  special 
properties  to  Sn+1  = X - Y that  are  the  source  of  the  rapid  convergence 
of  the  method. 

Figure  9 shows  a^block  requiring  a one-dimensional  minimization  of 
a*  of  the  function  4> (X  + a Sq) . The  one-dimensional  minimization  uses 
a four-point  cubic  interpolation.  It  finds  the  minimum  along  the 
direction  Sq,  where  X is  the  coordinate  of  the  previous  minimum.  By 
trial  and  error  it  finds  three  points  with  the  middle  one  less  than  the 
other  two.  It  makes  a quadratic  interpolation,  and  then  a cubic  inter- 
polation. If  the  actual  function  evaluated  at  the  new  interpolated 
point  is  not  sufficiently  close  to  that  of  the  preceding  point  or  if 
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The  objective  function  is  linearly  dependent  on  the  design  variables; 
however,  the  constraints  are  both  linearly  and  nonlinearly  related  to 
the  design  variables.  The  minimum  area  of  steel  is  a linear  constraint. 
Figures  11  and  12  show  the  shear  stress  and  the  deflection  as  being 
nonlinearly  related  to  the  thickness  of  the  concrete.  Note  that  the 
shear  stress  is  almost  linear  and  is  constant  (independent  of  thickness). 
Figure  13  shows  the  useable  region  bounded  by  flexure,  shear,  and 
minimum  steel  constraints.  The  optimum  least-cost  solution  is  shown. 

This  specific  example  solution  considers  an  unlaced  section;  thus,  the 
maximum  shear  constraint  is  active.  Laced  sections  eliminate  the  shear 
constraint.  If  the  number  of  sides  supported  were  increased  from  N=2 
to  N=3,  the  design  space  would  change  as  shown  in  Figure  14.  There  are 
two  regions  that  are  useable  areas.  Obviously,  the  lower  one  offers 
the  least  cost  and,  therefore,  is  more  desirable. 

There  is  clearly  a complex  interaction  of  constraints.  Unfortunate- 
ly, the  optimum  solution  found  by  the  program  depends  on  the  starting 
point  selected.  The  program  converges  on  the  closest  relative  optimum. 
Several  alternative  starting  points  should  be  used  to  verify  a question- 
able optimum.  Revising  the  design  parameters  could  possibly  shift  the 
constraints  such  that  only  one  useable  solution  would  appear.  However, 
a slight  increase  in  shear  stress  (10%)  can  significantly  reduce  cost 
by  allowing  the  near-optimum  nonfeasible  solution  to  be  accepted. 

The  dual-space  problem  of  finding  a useable  solution  is  limited  to 
unlaced  concrete  slabs  only  because  lacing  eliminates  the  shear  con- 
straint. Nonautomated  design  for  these  conditions  is  almost  impossible 
when  one  considers  the  complexity  of  the  design  space  and  the  large 
number  of  iterations  required  when  an  initial  solution  is  not  feasible. 

Cost  data  used  in  the  program  can  be  selected  by  the  user.  However, 
the  data  used  herein  is  based  on  work  by  Picatinny  Arsenal  on  contract 
with  Ammann  and  Whitney  [7].  Table  2 shows  a comparison  of  unlaced  and 
laced  concrete  walls  with  and  without  sand.  The  example  considers  a 
15-foot-high  by  12-foot-wide  wall  subjected  to  a 200-psi,  10-ms  tri- 
angular loading  function.  In  all  cases  the  laced  concrete  (12-degree 
rotation)  is  less  expensive  than  unlaced  (2-degree  rotation)  designs. 

The  costs  for  sand/concrete  composite  construction  are  for  only  the 
front  wall.  When  the  rear  wall  is  included,  the  costs  almost  double, 
thereby  making  this  form  of  construction  unsuitable  for  relatively 
low  pressure  loadings.  It  should  be  pointed  out  that,  for  the  N=3  and 
4 conditions,  the  optimum  design  selected  is  actually  a near  optimum 


with  the  shear  capacity  slightly  exceeded  as  shown  in  Figure  15. 
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Table  2.  Comparison  of  Optimum  Solutions 
(200  psi;  10  ms;  wall,  12L  x 15H) 


N 

Side 


Theta 

(°) 


Sand 

(in.) 


Cost 

($) 


N=2 


i 

>7777777777’ 


N=3 


'! t ////In 


N=4 

t til  Li 


ittt,tLty 


>11/ 1/71  •* 


2 

12 

2 

12 

2 

12 

2 

12 

2 

12 

2 

12 


0 

0 

24 

24 

0 

0 

24 

24 

0 

0 

24 

24 


3,290 

2,289 

2,209° 

1,856° 


2,753 

2,019 

1,944' 

1,943° 

2,001^ 
1,958 
2 , 001' 
1,943° 


a,2> 


^One  wall  only  in  composite  construction. 

■JL 

Shear  capacity  exceeded. 


The  program  contains  an  option  to  analyze  walls  with  openings. 
During  many  analyses  it  was  noted  that  blast  doors  with  resistances 
much  higher  than  those  of  the  walls  transfer  significant  reactions  to 
the  walls  such  that  the  walls  are  incapable  of  accepting  these  and  fail. 
Computational  problems  arise  in  the  program  when  this  happens  in  that 
yield  regions  cannot  be  brought  into  equilibrium  by  yield  analysis 
methods.  To  avoid  termination  of  the  solution  at  this  point  the  door 
resistance  is  reduced  automatically  by  a factor  of  2 to  reduce  the 
reaction.  This  usually  allows  for  a successful  termination.  Unfor- 
tunately, this  destroys  the  original  starting  point  for  optimization, 
and  creates  problems  when  a nonfeasible  low-cost  solution  is  lost  and 
cannot  be  used  to  provide  direction.  It  is,  therefore,  not  possible  to 
perform  optimization  solutions  of  walls  with  openings.  Generally,  it 
has  been  found  that  compatible  designs  occur  when  the  door  is  designed 
to  have  approximately  the  same  resistance  as  the  wall. 
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Figure  2b.  Wall  geometry 


Figure  4.  Computer  results,  Example 
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DYNAMIC  CONCRETE  STRENGTH  5000,00 

DYNAMIC  STEEL  STRESS  U8000.QO 

THICKNESS  CONCRETE  INCHES  2^.0000 

thickness  of  sand  inches  ae.ooou 

THETA  ALLOWABLE  DEGREES  12.0000 


AREA  VERT  TOP  STEEl/FT  I.S800  COVER  2,0000 

AREA  VERT  faOT  STEEL/FT  1.5800  COVER  2.0000 
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Figure  6.  Computer  results.  Example 


r 


•— 

z 

U 

r. 

»— 

3 


O O O O O ^ ^ D 
• O O O O m 4)  Al 
ooo  o o r- 

J)  • • • • «-4  *4  t 

p-4  ^4  • • flO 


H M « M N ••  U 

X X 

3 C 

^ or  *—  *- 

»-  O X u u 


z 

H- 

3 

< 

<4 

Or 

3 

Q 

•— 

r 

3 

1— 

X 

a. 

3 

•— 

-4 

3 

•a 

13 

f- 

3 

X 

V- 

a 

1) 

3 

3 

41 

3 

u 

X 

3 

3 

44 

X 

O K 

i 

< 

-Li 

4 

4» 

U- 

• X 

3 

►— 

X 

X 

•— 

3 

3 

3 

O i3 

_J 

o 

♦H 

% 

4* 

X 

CO 

CO 

UJ 

X 

iX»  M 

< 

O 

aJ 

>- 

r 

3 

UJ 

X UJ 

• 

* 

3 

3 

CO 

•-f 

3 

UJ 

44 

X 

* 

o 

X 

Jj 

X 

3 

3 

3 

J) 

aJ 

1- 

XI 

Ui 

z 

3 

Z 

CO 

3 

a 

<o 

It  »- 

IO 

j 

CO 

? 

•J 

co 

<4 

jU 

(Jj 

a 

■5 

X 

m 

Al 

in 

o 

A 

x 

3 

Al 

X 

3 

ill 

X 

ill 

CO 

h- 

UJ 

T 

X 

X 

a* 

AJ 

A* 

• 

• 

• 

• 

•O 

9 

4* 

3 

44 

cO 

ill 

CO 

X 

X 

X 

lO 

•u 

r* 

rO 

A) 

• 

.»  z 

• 

r 

o 

3 

X 

41 

X 

M 

•*4 

3 

X 

X 

a 

x» 

f\l 

.-4 

to 

r- 

a 

X 

X 

-1  Cl 

o 

u 

4* 

I 

13 

*3  X 

3 

H— 

z 

X 

ill 

N_a 

X 

Al 

0 

• 

w M 

o 

3 

■> 

.3 

►—  H- 

XI 

• 

z 

<3 

X *-i 

Z 

no 

o 

Al 

CJD  CO 

«-* 

AJ 

-•  3 

• 

ill 

•— 

co 

UJ 

* E 

It 

2 

c/> 

o 

AJ 

4-X 

o cr 

XI 

Xi 

o 

o 

or 

o 

3 

UJ  U 

r 

A* 

o 

o 

a 

O 

o • 

o 

X 

• 

• 

• 

• 

•-4 

AJ 

vn  »o 

3 

o 

w 

• 

x 

o 

3 tO 

M 

X 

CO 

X 

X 

• 

• 

• 

• 

X 

a. 

ar  iii 

• 

X 

o 

•-4 

• 

• 

• ^ 

3 

X 

a 

9 

— ♦ 

X 

AJ 

to 

X 

• 

4 > 

x 

iTI 

o 

X f\i 

z 

UJ 

X 

•>4 

X 

in 

-o 

t* 

• 

x -< 

<4 

X 

• 

• 

rvi 

• 

•— * 

>• 

O 

u co 

o 

o 

O 

o 

3 

• 3 

o 

r* 

1 

1 

f 

• -1 

1*0 

m n 

X 

■ a 

• 

N 

M 

II 

M 

M 

it 

H 

X 

• X 

z 

XI 

X 

111 

• UJ 

3 

<4 

•— 

X 

D 

o 

HH 

3 

M 

*s 

X 

< 

co 

J\  X 

fO 

o 

Al  ^ 

aJ  x 

<3 

o 

—4 

cO 

►—  X 

44 

3 

Ai 

in 

X 

X 

o 

*4 

c 

*■4 

% 

*-♦  cr 

x. 

3T 

a 

<X 

X 

X 3 

XI 

3 

o 

XI 

o 

o 

aj  po  m 

ill 

*-  3 

3 

X ) 

X 

-1 

41 

>-> 

XI 

3 

CO 

O 

X 

X 

• 

• 

• 

Jr  a. 

<x 

• 

3 

co 

3 

X 

X 

ill 

XI 

T 

o 

X 

• 

• 

• 

• 

• 

• 

•4 

*4 

•4 

CO 

m cil 

u 

c 

U 

TC 

-1 

UJ 

X 

3 

►-# 

X 

o 

a 

• 

-1 

aJ 

%mS 

■X 

*- 

3 

o 

*-> 

7 

co 

< 

3 

i 

<4 

t— 

t— 

i 

CO 

— * 

1- 

44 

at 

o 

UJ 

T 

3 

CO 

X 

X 

►- 

3 

Ol  *- 

o 

H 

y. 

j> 

Z 

'll 

ill 

X 

X 

tO 

X 

< 

o 

*-% 

ill 

< 

M 

o 

X 

UJ 

w 

*-4 

Xl 

z 

aJ  3 

o 

4l 

I> 

r 

111 

aJ 

G 

cr 

a 

H 

a 

p— 

3 

•— 

> UJ 

• 

M 

<a 

4 

3L 

:> 

M 

< 

Vli 

X 

M 

X 

o 

X 

o 

--4 

»o 

in  n 

M 

«-4 

jy 

* 

X 

•-« 

H— 

i 

X 

X 

3 

o 

■*4 

CO 

3 

rvi  iO 

AJ  X 

fO 

X 

to  X 

Kl 

i) 

CO  JL 

a 

111 

X 

<* 

3 

UJ 

UJ 

3 

UJ 

o 

X 

Ai  rvi 

•-4 

AJ 

9 

in 

X 

o 

3 tal 

Sm* 

* 

t- 

3 

3 

X 

v 

X 

x> 

1- 

X 

XI 

3 

CO 

in 

X 

X 

• 

• 

• 

• 

• 

9 

• 

«li 

3 <E 

•-4 

X) 

o 

3 

X 

3 

XI 

r 

X 

M 

M 

X 

a 

X 

• 

• 

9 

«*4 

^4 

*-• 

•4 

AJ 

a x 

T 

3 

X 

4 

a 

Q 

cO 

44 

44 

►- 

CO 

H 

X 

H- 

z 

X 3 

< 

X 

z 

X 

X 

V, 

< 

X 

X 

j 

ill 

1“ 

XI 

aJ  Z 

X 

in 

HI 

u 

111 

3 

u u 

3 

44 

a 

iwpuLse  (psi.^stc) 


o c*  o 

o o c 


o o o 

O 3 lT* 


i o fvj  ,r  os  h 

o 

o 


ju  »i!  « » ^ ^ ^ ^ 'Vi  m o n 
i a:  m a 3 <vi  n»  • • • • *n 

e x a.  x»  cm  -^  ~*  **>  a t>  x> 

X*a>w  S ^ f\!  ^ O 

o > 

Z o 


2 W 

UJ  UJ  O o JV  ^ S n 3 ^ fv  o 

ox*-*  • • • • « 'M  ^ in  »n 

m a.  in  ^ 3 ® ® • • • • x> 

u (t  a • 

z ujwapn^'**»n’n^N  • 
•-«  > o 

3 


u ai^inxxon^o^ 
* uj  Dino.noo?\jfOin 

u u)  ^ >o  s*  ^ •••• 
a r ••••••—  — — — 

t ^ o 
(O 


il  m o s>  o ® o ^ *n  3 in  n 

4ioumiPs<\i®rn«fnoM 

Qii)'xinih*^(\i3inN®o 

ni  j </)  in  x ••••••• 

i a.  x.  • • • —*•-*  »-•  m 

M x ^ 

s-  ui 


K X S 
lL  T IL 


■ — 4) 

i o m 


CO 

XNII 
• O 

»-*  >-  UJ 
« 2 *- 
a ui  o 

W -D  4jJ 

M 3 
UiUlL 
05  2 ui 
-|m  a 
3 
<L 
£ 


ui 

C5  s— 
X X 
4 J 
X — * 

U K Hi 
X X 
tL  19 
O — j 
aJ  3 
ui  * < 
U I 
2 Ui 
•*(!>►- 
k a ® 

« 4 4 
M I J 

ou® 


CO 

o 

• • • 

X 

Ui 

4> 

*■“  *•  ^ 

• 

X 

•— 

liU.lL 

M 

X 

• 

CO 

fO 

X 

in 

*\i 

o 

«-4 

(X 

4) 

3 

• 

«-• 

3 

<\i 

o 

• 

* 

»-• 

IT. 

r*> 

UJ 

• 

X 

3 

o 

(O 

G <* 

CO 

2 

Ui 

3 

X 

3 

CL 

X UJ 

o 

O C5 

< 

X 

a 

Ui 

3 

< 

JJ  <1 

CO 

3 

Ui 

> X 

3 

a 

3 O 

3 

U 

cc 

a 

3 

03 

4 Z 

X 

3 

X Ui 

M 

Z 

O 

*-  UJ  UJ 

O 

*— 4 

'JC3  1 

Ui 

*-• 

*— 

Z X l- 

c 

3 

»— 

•—4 

lil  4 ill 

3 

-4 

< 

L- 

3 X CD 

u 

►- 

ar 

o 

O 

a 

3 

*-4 

3 • 

2 

*- 

Q 

u 

3 u.  *-  O 

4 O ® M 
2 *-•  ►- 

»-  o o 
•-  x U1 

03  LD  • 3 

4 M Z li. 
J UJ  M UJ 

® I r or 


DYNAMIC  CLiNCHfcTE  STRENGTH  5000,00 

Dynamic  STEEL  STRESS  48000. 0U 


o o o o 
o o o o 
0 0 0 3 
0 3 0 0 


a*  <\i  m t+i 


X X x <x 


\ N 

£j  L •/)  0) 

j ji  a. 


Ui 

UI  Ui  UJ 

t 

> > > 

o 

> 

CO 

CO 

•-I 

3 3 3 3 

Ob 

a 

X 

3 

*n 

u 

u uu 

• 

• 

• 

• 

fM  a- 

—4 

X 

© 

XI 

rw 

4) 

3 

o 

3 

3 

X 

OJ 

X 

© 

3 

3 3 3 

3 

3 0 3 

« © 

3 

O i<> 

«•* 

O 

3 

o 

3 

3 

X 

3 

Al 

9 

1 

0 3 0 

3 

O O O 

a 

ON" 

4) 

• 

• 

ui 

o o 

m 

• 

• 

• iT 

in  An 

0 3 0 

*> 

X)  JD  X) 

• 

• 

■ • 

O 

in 

o 

a 

o 

3 

3 

O 

O 

3 0 0 

jun  ioit 

f\l 

Nl  l\l  o 

•-4 

l\i 

• 

• 

© 

© 

T 

MO 

• • • 

• 

• • • 

in  r*i  ^ 

• 

3 

3 

• 

^ 3 M 

^ 

*-# 

*)  r-* 

o 

3 

ao 

AO 

o 

►— 

tf\l 

4-* 

ul 

X 

X 

w+ 

4) 

X 

-1 

O 

w— 

«-4 

Zt  *- ; 

”M  f\J 

X 

X 

x> 

■u 

w— 

w$ 

X 

o- 

il 

*— 

A— 

zj  a. 

in  <o  40 

X 

Ui 

A— 

A- 

JL 

X 

►— 

3 

X 

A- 

73 

& 

ui  uJ  aU 

X. 

z 

2 

3 

►- 

r 

X 

T)  X 

3 

L I Ui 

XI 

u 

A- 

A— 

A~ 

ui 

Ul 

3 

-m 

3 

X 

3 

X 

(O 

UU  1 

r—  A— 

m; 

ul 

2 

2 

X 

X 

3 

2 

X 

X 

3 

z z e 

1- 

•“  llL  L. 

UJ 

•s 

Li 

A— 

u UJ 

3 

u 

u 

Ui 

> 

3 

X 

X 

X 

*-#  — a ui 

U. 

U.  N X 

X 

►— 

3 

X 

X 

ul 

T 

X 

X 

X 

J T 

*— 

o 

3 

X 

3 

Q 

X. 

V -l-l 

— • u 

-a 

X 

Ui 

3 

3 >- 

3 

o 

UJ 

•—A 

to 

X 

-1 

J Ui  UJ 

m 2 

M 

uJ 

ui 

(L 

a 

Ui 

r 

T 

-J 

3 

CO 

> 

Ui 

ui 

3 

A— 

It. 

ui 

ui 

ui  Lti  UJ 

a 3 

Z 4 

a— 

Z 

a 

X 

<1 

Ui 

3 

2 

•a 

< 

A— 

7) 

•—  a 

UJ 

Uir-H 

u 

X 

M »H  tfj 

fO 

•—a 

3 

LI 

-j 

3 *“ 

3 

© 

»— • 

X 

<1 

M 

ui  z 

K 0)  U) 

in  \ 

ui 

A— 

• 

ui 

<o 

z 

< 

< 

2 

2 

•— • 

4 

3 

3 

4 

3 

o <o  1 

IX  •*.  ui 

0»  tf) 

3 -1 

Z 

X X 

a— 

© 

3 

X 

3 

LI 

Li 

3 

3 

<© 

3 

< 

a 

— 

UlO  J 

a h 

3 Ui 

•— • 

O ul 

z 

V. 

z 

>— 

LI 

M 

N 

O 

Cj 

•3 

•j 

A— 

a 

z a 

a i-  0 3 

3 UI 

X 2 

ui 

N 

3 

«* 

A— 

A- 

M 

H 

**> 

>- 

3 

3 

O 

_» a 

•a 

3 U.  •«* 

3 0 1-® 

O »— 

»— 

M 

u 

M 

3 

3 

X 

X 

X 

or 

Ui 

Ui 

u 

3 

►— 

UO  J 

XI 

3 » 

X 

3 

X 

3 

X 

X 

X 

UJ 

Ui 

Ui 

3 

3 

»-• 

M 

3 Z 

X 

< 

3 

•M  rs* 

XI 

Ui 

UI  A- 

ul 

T. 

■a 

<L 

CJ 

■» 

> 

r 

r 

2 

Ui 

>■ 

► 

3 

<* 

X 

< 

tr 

ifl  n j 

r— 

K *-<  M 

/> 

r 

X 2 

X 

ul 

ill  TC 

3 

2 

Ui 

ui 

< 

tf)  lO  J 

X 

x x a 

ui  3 

Cl  U 111 

ac 

X 

X 

X 

3 

Ui  UJ 

ul 

ui 

M 

Z 

2 

<3 

X 

ui 

r 

Ul 

UJ  UJ  41 

jJaiOC 

f-  X 

c 

< X 

-D  3 D 

in 

J)  UL- 

■-» 

-> 

> 

» 

L- 

3 

3 

3 

A— 

n 

r 

CO 

c 

2 Z 

> 

>•  X X 

Ui 

X 3 

(J 

A— 

z 

A** 

M 

M 

>-# 

X 

M 

»-• 

CO 

in 

< < 4 

X 3 

to  «J  X 

o u 

X 

X 

A— 

A— 

A— 

A— 

3 

3 

A— 

A~ 

X 

MnJ 

MM 

OUK 

«• 

< u < 

U M CO 

X 

< 

3 

3 UJ 

M 

•u 

M 

<1 

X 

3 

< 

< 

A-A 

*— * 

1- 

*— A 

*- 

UHllI 

ill  Ui  Ui  Ui 

z *- 

D UJ  ai  HI  Ll 

U. 

-1 

J X 

CO 

3 

to 

<3 

X 

ui 

L> 

3 

A- 

X 

X 

X 

3C 

XXX 

X 

XXX 

a < 

X 

> > 

> 

-J 

-I 

2 

3 

ui 

3 

ui 

3 

A-# 

r-. 

3 

3 O Ui 

3 

UJ 

*< 

«*•  < ^ 

U X 

3 

C *< 

<* 

a a 

•a 

<* 

3 

X 

Z 

X 

2 

*© 

>“ 

a _»  a 

X 

•> 

X 

X X 


o o 

»1  M 

X z 

z z 


C D i/)  (O 

j a a 


-d  a*  »/>  'xj 


^ ^ x» 


f 

o 

o 

o 

*M 

x> 

3 

o 

i 

o 

3 

3 

t 

X 

>1 

a 

■4 

'V 

Ol 

X iO 

a 

o 

CM 

%% 

• 

• 

m 

• 

IT 

in 

4 

4 ao 

• 

r-  o 

3 

O 

O 

^“4 

O • 

-4 

a 

o ^ 

ai 

, • • 

o 

o 

t> 

r*> 

XT 

A CD 

o 

• «4 

• CM 

f o o 

r* 

• 

• 

• 4> 

• 

o • 

o 

• 

o o 

•o 

4> 

X> 

•— 

A 

•“U 

>o 

X 

*4 

<T  o 

3 

6— 

ru  fM 

X 

X 

— • — • 

a 

3 

»— 

3 

a 

»- 

►— 

JL 

X 

A— 

3 

a 

»- 

03 

a 

2 

a 

3 

*- 

X 

X 

O Jt 

3 

z 

K *- 

uj 

Ui 

D 

U> 

X 

3 

3L 

A 

z 

o 

z z 

X 

£ 

3 

Z 

•4 

X 

X 

3 

3 

A— 

•4 

Itlui  3 

-w 

a. 

Ur 

>- 

3 

X 

X 

X 

•4 

a 

UJ 

X X 

X 

X 

3 

X 

— 7 

u. 

3 

♦— 

CO 

X 

3 

3 3 

uj 

•4 

tf> 

•X 

U 

A— 

a 

UI 

X X 

3 

3 

IO 

> 

vU 

uJ 

o »- 

►- 

u. 

U J 

a 

4 

*4  3 

VtX 

4 

4 

UJ 

a 

z 

z 

4 

4 

•— 

J> 

3 

Ui 

>- 

ui 

z 

X U 

u 

Z 

-1  3 

w 

3 

o 

*— * 

*4 

X 

4 

M 

►— 

u. 

•3 

3 

»— 

X 

3 

a 

•4  Ui 

z 

3 

* 4 

z 

a 

«4 

4 

3 

-I 

4 

O 

o 

A 

uJ 

ui 

X z 

z 

•4 

3 O 

4 

4 

u u 

3 

3 

CO 

u 

4 

o 

M 

O 

ui 

3 3 

M 

»— 

f- 

M M 

M *4 

O 

o 

O 

<4 

►- 

e 

3 

►- 

3 

»— 

o 

4 

O O 

CO 

o 

K 

M 

*4 

>- 

_J 

3 

o 

-i  a 

4 

X 

X 

Cl  ui 

a 

•4 

ui 

64  4 

•4 

Ui 

or  at 

•x 

'V 

UJ 

4 

3 

*— 

4 

Ui 

4 f— 

X 

3 

O H K 

CO 

3 

UJ  Ui 

3 

■mj 

*4 

M 

3 

X 

JL 

4 

X 

U.  UI 

3 

*4 

u. 

3 A A 

Ui 

u. 

■>  > 

X 

t 

z 

<4l 

v- 

V 

J 4 

X 

4 

X 

>- 

X 

3 

_l  Ui 

Ui  4 4 

a 

Ui 

3 

z 

UU 

4 

Ui 

4 

3 

*— 

CO  o 

ui 

Q 

•4 

_j  _J 

a 

UI  UJ  1 ii  UJ 

M 

z 

z 

tu 

X 

Ul 

X 

Ui 

3 

H 

CO  z 

M 

>- 

cl  a 

UJ 

'>■  > 

> 

> 

3 

3 

3 

#» 

(A 

X 

CO 

X 

O 

U 

4 3 

>- 

3 

3 

< 

i— 

o 

P-*  M 

*4 

or 

M 

M 

4 

CO 

A 

4 

4 

X u 

Hi 

•4 

o o o 

4 

♦4 

6—  6- 

6— 

K- 

3 

3 

»— 

►— 

X 

>sj 

<NI 

* 

a 

♦— 

*— 

z 

K K 

X 

►- 

M 4 

M 

4 

X 

3 

4 

4 

1-4 

M 

1- 

*4 

»- 

3 

4 

o co 

CO 

CO 

(O 

O 40  40 

4 

co 

CO  *3 

CO 

f3 

X 

Ui 

u 

u 

►- 

X 

X 

r 

X 

3 

U 

4 CO 

X 

4 

4 

0 4 4 

►“ 

4 

O ui 

3 

Ul 

3 

M 

c 

3 

3 

a 

Ui 

3 

Ui 

3 

C 4 

M 

3 

3 

UJ  _l  _l  _l 

3 

X z 

X 

a 

yf> 

>- 

-1 

-1 

3 

X 

■> 

X 

.> 

4 

X 

3 X 

u. 

ui 

ui 

to 

JJ  ui 

a 

aJ 

X 

CO 


26 


ULTIMATE  RESISTANCE  RU  66. a* 

ELASTIC  DEFLECTION  LIMIT  XE  .1775 

STIFFNESS  KE  37a, 50 


i 


BOOOOOOOOOOOOOOOOOOOOOOOOO 

&oo  ooooooooooooooooooooooo 

9000C000030000000000000000 

Boo  ooooooooooooooooooooooo 

BOOOOOOOOOOOOOOOOOOOOOOOOO 


C K 

r- 

x> 

9 

in 

^4 

in 

O' 

r\i 

at 

in 

X 

X 

-0 

a 

Oi 

O' 

X 

f\i 

- 

X 

o 

K1 

rw 

x 

« 

9 

X 

•-• 

O' 

91 

© 

■n 

*V 

o 

4\l 

m 

o rv 

m 

x 

p^ 

4) 

X 

x> 

m 

'■n 

p- 

*\l 

X 

a 

rv 

X 

r- 

X 

X 

9 tf> 

X 

X* 

a 

o 

r\j 

rO 

.n 

x 

x 

X 

X 

D 

9 

a 

9 

a 

9 

O' 

p r- 

p^- 

p^ 

«o 

CO  CD 

X f 

X 

X 

X 

X X 

X 

X 

X 

X 

X 

X 

cac 

« to 

9-4 

•** 

*-• 

9-$ 

*4 

•-• 

«-• 

•-* 

9-4 

9-4 

4-+ 

•-• 

— • 

9-4 

9+ 

MVK<vi'ryc^oxu)'<)<tr>iff  3 a 3 9 ^oinciroi/i 
^3oxx'\ini'<ixi33irin<c^r>-<ic90'Oxx(M»i 
nnAi'<o»<'r>-«i<ta  rt  ft)  •«  o <mc  r-  c lP  3 3 ift  fti  - > o 
ttftlftl'Mftl'-'*^'*''— '•’-<'<-•3  030030000  O 


BOSflC9SfCllC(IC9X>rit1ITirflKitlSfC 

kr.r-r-r-r<Kr<r-NNNr’r<KKr<r<r<NKNNNNr' 

B a O 30000000000000  0000  00000 


9 <S 

► 3 

.*  p- 


9 9 


-^iPs-^irsftiirxiftiioffftiirvfti  cffivooiftOB 
O 1/1  afi-OftlNftlti'ftS  99  3 O /I  O a -I  « ft)  N ft) 
3BlTO)i—r’ft)«J'ftB3  0-4>«-*^-<\)30>*»99  o 
— 3®-<irBftlB»iftroiftr'-<3  X>  — I IT  ® ft)  IP  (ft 
CC  — • 3 D — • 3 X -*3®-’l)lC-<ir®— 

3 O^OBmBIVMVCKU)  3 9 3 0 1/10  O - ftl  K 


•o 

X 


D 

r 


1/1  O o N N X>  D 9 9 O O -<  -•  ft)  ft)  1ft  3 3 1/1/1®  D 
l»tlftl»t>ftift«llftFftKI333333  3333333 


m 

a 

X 

o 

X 

X 

o 

rv 

in 

X 

o 

• 

—4 

© 

X 

*n 

in 

-n 

O X 

'v 

n 

mi 

rv 

X 

N ^ 

m 

X 

•n 

J> 

in 

• • 

X 

a 

rv 

X 

u 

O X 

• 

• 

• 

• 

• 

• 

ul 

X 

X 

X 

«4 

9-* 

X 

»“ 

14. 

>» 

z 

P— 

o 

M 

#— • 

p- 

u 

a 

•—i 

3 

o 

3 

X 

3 

ul 

M 

•— » 

Ul 

>- 

-1 

X 

3 

> 

p— 

u. 

Ui 

2. 

ul 

a 

2 

3 

3 

o 

3 

3 

3 ► 

$-* 

3 

*-* 

■V 

3 P- 

X 

X 

Ul 

p- 

X 

■ 1 9-4 

c.- 

a 

3 

X 

u 

u 

X 

tf>  > u 

3 

ui 

X 

3 

z 

Ul 

X 3 

»— 4 

3 

»-* 

P- 

X 

3 

t- 

X X 3 

X 

14. 

X 

<. 

P— 

iL 

2 

a < jj 

u 

UI 

< 

2 

X 

Ui 

ul 

4 u > 

a. 

3 

X 

X. 

*— • 

O 

X 

3 X 

2 

X 

;3 

3 X 

-i 

X 

3 

a 

Ul 

u 

X 

C uJ  X 

3f 

3 

P- 

M 

X 

M 

1 

OJ  u 

X 

X 

P— 

*v 

p- 

14. 

< X 

3 

*—♦ 

ui 

<x 

o 

X 

3 or 

P— 

x 

X 

X 

< 

•4 

X 

3 Ul  X 

< 

•v 

*-* 

3 

u 

31 

< 

^ < 

2 

X 

Q 

3 

Ul 

X 

1 4 X 

Start 


Figure  9.  Logic  diagram  for  minimization  of  (}> (X) 
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Figure  15.  Revised  design  space,  N=3. 


Appendix 


EXAMPLE  PROBLEMS 

This  Appendix  presents  several  example  problems  showing  the  use 
of  optimization.  Test  Case  A-l  is  shown  in  Figures  A-l  and  A-2.  The 
input  estimated  24-inch-thick  concrete  and  1.58  sq  in.  of  steel  per 
foot  with  a cost  of  $8,199.  The  final  design  resulted  in  a 12- inch-thick 
concrete  section  with  vertical  steel  0.41  sq  in. /ft  and  horizontal 
steel  1.23  sq  in. /ft  costing  $4,958. 

Input  for  Test  Cast  A-2  is  shown  in  Figure  A-3.  The  preliminary 
cost  was  $3,433;  however,  the  shear  capacity  of  the  proposed  section 
was  inadequate.  The  final  design  was  a 23.6-inch-thick  concrete  section 
with  0.86  sq  in. /ft  of  vertical  steel  and  0.98  sq  in. /ft  of  horizontal 
steel  costing  $3,192,  as  shown  in  Figure  A-4. 
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Figure  A-l.  Input  for  Test  Case  A-l. 


Figure  A-2.  Results  for  Test  Case 
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Figure  A-2.  Continued 
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Figure  A-3.  Input  for  Test  Case 
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Figure  A-4.  Continued 
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NAVFACENGCOM  • SOUTH  DIV.  Code  90.  RDTAEl.O.  Charleston  SC;  Dir..  New  Orleans  LA 
NAVFACENGCOM  - WES  I DIV.  408,  San  Bruno  CA;  AROICC.  Point  Mugu  ( A;  Code  (MB;  Codes  09PA 
NAVFACENGCOM  CONTRACT  Bethesda.  Design  Div.  (R.  Lowe)  Alexandria  VA:  Eng  Div  dir.  Southwest  Pac. 
Manila.  PI;  OICC/ROICC.  Balboa  Canal  Zone;  ROICC  l.ANT  DIV..  Norfolk  VA;  ROICC.  Pacific.  San  Bruno  CA; 
TRIDENT  (CDR  J.R.  Jacobsen).  Bremerton  WA  98310 
NAVMARCORF.STRANCEN  ORU  1 1 18  (Cdr  D.R.  Lawson).  Denver  CO 
NAVMIRO  OIC.  Philadelphia  PA 

NAVOCEANSYSCEN  Code  6565  (Tech.  Lib  ).  San  Diego  CA 
NAVORDSTA  PWO.  Louisville  KY 

NAVPHIBASF.  CO.  AC'B  2 Norfolk.  VA;  UCT  I (MacDougal)  Norfolk.  VA 
NAVREGMEDCEN  SCE  (LCDR  B.  E.  Thurston).  San  Diego  CA 

NAVSCOLCECOFF  C44A  (R.  Chittenden).  Port  Hueneme  CA;  CO.  Code  C44A  Port  Hueneme.  CA 
NAVSEC  Code  6034  (Library).  Washington  IX' 

NAVSECGRU  ACT  PWO.  Torri  Sta.  Okinawa 
NAVSHIPREPFAC  Library.  Guam 

NAVSHIPYD  Code  202.4.  Long  Beach  C'A;  Code  400.  Puget  Sound;  Code  410.  Mare  Is..  Vallejo  CA;  Code  440 

Portsmouth  NH;  Code  440.  Norfolk;  Code  440.4.  Charleston  SC;  Library.  Portsmouth  NH;  PWDll.T  N.B.  Hall). 
Long  Beach  CA:  PWO  Portsmouth.  NH;  PWO.  Mare  Is. 

NAVSTA  CO  Naval  Station.  Mayport  FL;  CO  Roosevelt  Roads  P.R.;  Engr.  Dir..  Rota  Spain;  Maint.  Div.  Dir/Code 
531.  Rodman  Canal  Zone;  PWD/Engr.  Div.  Puerto  Rico;  PWO.  Keflavik  Iceland:  PWO.  Mayport  FL;  PWO.  Puerto 
Rico:  SCE.  Guam;  SCE.  Subic  Bay.  R.P.:  Utilities  Engr  Off.  (LTJG  A S.  Ritchie).  Rota  Spain 
NAVSUPPAC'T  AROICC  (LT  R.G.  Hocker).  Naples  Italy;  CO.  Seattle  WA;  Engr.  Div.  (F.  Moll  tea).  Naples  Italy 
NAVSUPSYSCOM  Code  NSUP-0323.  Washington  IX 
NAVSURFWPNC'F.N  PWO.  While  Oak.  Silver  Spring.  MD 
NAVTECH  I RAC'EN  SCE.  Pensacola  FL 
NAVTRAEQUIPCKN  Technical  Library.  Orlando  Fl. 

NAVWPNCEN  PWO  (Code  70).  China  Lake  CA;  ROICC  (Code  702).  China  Lake  CA 
\ A VWPNENGSU  PPACT  Code  ESA-1163.  Washington  IX' 

NAVWPNSTA  ENS  G.A.  Lowry.  Fallhrook  C'A;  PW  Office  (Code  090  ) Yorktown.  VA:  PWO  Yorktown.  VA 
NAVWPNSUPPC'EN  Code  09  (Boennighausen)  Crane  IN 
PM  1C  Code  4253-3.  Point  Mugu  CA 
NAVEDTRAPRODEVCEN  Tech.  Library 

NAVFACENGCOM  - l.ANT  DIV.  Eur.  BR  Deputy  Dir.  Naples  Italy 
NAVSUBASE  LTJG  D.W.  Peck.  Groton.  CT 

NCBC  CEL  (CAPT  N.  W.  Petersen).  Port  Hueneme.  CA;  CEL  AOIC  Port  Hueneme  CA;  Code  10  Davisvifle.  Kf:  PW 
Engrg.  Gulfport  MS:  PWO  (Code  80)  Port  Hueneme.  CA 
NCR  20  Code  R3I;  20.  Commander 
NMCB  5.  Operations  Dept.;  THREE.  Operations  Off. 

NRI.  Code  84(8)  (J.  Walsh).  Washington  DC:  Code  8441  (R.A.  Skop).  Washington  IX' 

NSC  Code  54.1  (Wynne).  Norfolk  VA 
N IC  Commander  Orlando.  FL 

NUSC  Code  EAI23  (R.S.  Munn).  New  London  CT;  Code  TAI3I  (G.  De  la  Cruz).  New  London  CT 
ONR  Code  7(8)F  Arlington  VA;  Dr.  A.  Laufer.  Pasadena  CA 
PMTC  Pat.  Counsel.  Point  Mugu  C'A 

PWC  ACE  Office  (LTJG  St.  Germain)  Norfolk  VA;  CO.  Great  Lakes  II  ; Code  120.  Oakland  CA;  Code  I20C  (A. 
Adams)  San  Diego.  CA:  Code  I20C  (Library)  San  Diego.  CA;  Code  2(8).  (ireat  Lakes  II.;  Code  220.1.  Norfolk  VA: 
Code  30C  (Boettcher)  San  Diego.  CA;  OIC  CBU-405.  San  Diego  CA:  XO  Oakland.  CA 
USAF  SCHOOL  OF  AEROSPACE  MEDICINE  Hyperbaric  Medicine  Div.  Brooks  AFB.  TX 
USCO  (G-ECV/61 ) (Burkhart)  Washington.  IX';  G-EOE-4/61  (T.  iX>wd).  Washington  IX:  MMT-4.  Washington  IX' 
USCC  RAD  C E NTER  Tech  Dir  Groton.  CT 

l SNA  Ch.  Mech  E ngr.  Dept  Annapolis  MD:  PWD  Engr.  Div  <C.  Bradfordl  Annapolis  MD:  PWO  Annapolis  MD 
LE  HIGH  UNIVERSITY  Bethlehem  PA  (Frit/  Engr  Lab  No.  13.  Beedle) 

LIBRARY  OF  CONGRESS  W ASHINGTON.  IX  (SCIENCES  A TECH  DIV) 

MASSAC  HUSF  I IS  INST.  OF  I FCHNOI.OGY  Cambridge  MA  (Rm  10-5(8).  lech  Reports.  Engr.  Lib  ):  Cambridge 
MA(Rm  14  E2I0.  Tech  Report  Lib.):  Cambridge  MA  (Whitman) 

SOUTHW  EST  RSCH  I N S I R DeHart.  San  Antonio  TX 

l S MERC  HAN  I MARINI  ACADEMY  KINGS  POINT.  NY  (REPRINT  CUSTODI  AN  I 
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UNIVERSITY  OF  ll.I.INOIS  URBANA.  11.  (LIBRARY) 

UNIVERSITY  OF  TEX  AS  Inst.  Marine  Sci  (l  ibrary).  Pori  Aransas  TX 
CONCRETE  TECHNOLOGY  CORP.  TACOMA.  WA  (ANDERSON) 

NEWPORT  NEWS  SHIPBI  IHi  & DRYIXX  KCO.  Newport  News  VA  ( lech.  I.ib.) 

NORWAY  Norwegian  lech  Untv  (Brandl/ueg).  Trondheim 

PACIFIC  MARINE  TECHNOLOGY  LONG  BEACH.  CA  (WAGNER) 


